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Abstract. We investigate the ultradiscrete KdV equation with periodic boundary con- 
ditions where the two parameters (capacity of the boxes and that of the carrier) are 
arbitrary integers. We give a criterion to allow a periodic boundary condition when 
initial states take arbitrary integer values. Conserved quantities are constructed for the 
periodic systems. Construction of background solutions of the periodic ultradiscrete 
KdV equation from the Jacobi theta function is also presented. 
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^ ■ 1. Preface 

^ ! In this paper we investigate the boundary conditions and the background solutions of 

I a cellular automaton which is called the Box Ball System (BBS). 

The content of this paper is as follows. The BBS is derived from the discrete KdV 
equation by a limiting procedure called 'ultradiscretization', which we will explain briefly 
in section 2. In section 3, we give a criterion for the BBS and the BBS with a Carrier 
(BBSC) to allow a periodic boundary condition even when the system allows 'negative 
solitons'. We also define the background solution of the BBSC using the conserved quan- 
tities of the BBSC. In section 4, we introduce a gauge transformed discrete KdV equation 
CN ! whose solutions converge to the upward-shifted solutions of the BBSC by ultradiscretiza- 
tion. We then investigate the relation between the Jacobi theta solutions of the discrete 



■ KdV equation and the background solutions of the BBS in section 5. In section 6, we 
. present an example of multi-soliton solutions with the boundary condition in which the 
! value in n — )■ oo and that in n — )> — oo differs from each other. 

O 

2. Discrete KdV equation and BBS 
2.1. Ultradiscretization. The discrete KdV equation is defined as 

where n and t take only integer values. The discrete KdV equation is transformed to the 
bilinear form 

(2.2) (1 + 5)a*+\a^i = 5at,\a'^^ + a* 



by putting 

w. 



t '^n^n+l 



The ultradiscretization is a limiting procedure in which the dependent variables of the 
discrete equations become also discretized [I]. The ultradiscretization transforms discrete 
equations into piecewise linear equations. First we use the following lemma to ultradis- 
cretize fl21|l . 
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Lemma 2.1. Under the boundary condition lim„_=._oo = 1, the discrete KdV equation 
(12. ip is turned into 



(2.3) 



w. 



t+i 

n+l 



w 



t+i 



Then we put w^^ 
obtain 



exp 



^) , 6 = exp (-j) {e > 0) and take the hmit e +0 to 



(2.4) U'^Xl = min L - f/*^,, i^l - U-^ 



k 



under the boundary condition hm„_)._oo = 0. 

We define the ultradiscrete KdV equation by the equation (12.41) ^2]. It is also equivalent 
to the time evolution equation of BBS with box capacity L [3, 4]. Next we ultradiscretize 
(12. 2p . Putting (T* = e'^"''^, 6 = e~^l^ and taking the limit e — > +0, we obtain the bilinear 
form of the ultradiscrete KdV equation: 

(2.5) r*+l + = max[r*;l + r*+i - L, r* + r*+i]. 

The ultradiscretization preserves the solitonic nature of the continuous KdV equation. 

Example 2.1. Let L = 1, then equation (12. 4p is closed under f/* G {0, 1}. We give an 
example of the time evolution below (where a dot indicates a zero). 



t=l 
t=2 
t=3 
t=4 
t=5 



111. . .11. . .1 

. . .111. .11. .1 

11. .11.11 

11. .1. .111 

11.1 111. 



2.2. BBS with a carrier. The time evolution of the BBS with a Carrier (BBSC) is 
expressed as follows [5J. We prepare the "carrier" which can carry at most / balls. At 
each time step in the evolution, the carrier moves from the left to the right. While 
the carrier passes the j-th box, the following action occurs. Assume that the carrier 
carries c (0 < c < /) balls before it passes the j-th box, and also assume that there are 
U {0 < U < L) balls in the j-th box. Then, when the carrier passes the box, the carrier 
puts min(c, L — U) balls into the box and receives min(f/, I — c) balls from the box. 
That is to say, the carrier puts as many balls into the box as possible and simultaneously 
obtain as many balls from the box as possible. This rule can be expressed in the following 
formula 



n— 1 \ / n n—1 

t+1 
k 



(2.6) f/*+i=min L-[/*, ([/* - ^*+i) + max 0, ^ C/* - 

\ k=—oo / \ k=—oo k=— 

under the boundary condition 

lim Ut = 0. 
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3. Periodic BBSC that allows negative solitons 

3.1. Negative solitons. Negative solitons are non-solitonic trains of negative values at 
a speed of 1. These arise for example when we put negative initial values to the BBS(C) 
equation. 

Example 3.1. An example of a negative soliton in a BBS with box capacity 1 is as follows 
(where a dot indicates a zero). 

t=l: . .-1-1-1 . .-1 .... 1 1 

t=2: . . .-1-1-1 . .-1 11... 

t=3: . . . .-1-1-1 . .-1 11. 

While the sequence of Vs has a speed equal to its length, the sequences of -1 have speed 1 
regardless of their lengths. □ 

The pioneering works on negative solitons and background solutions are due to Hirota 
[6], and Willox et.al. [7]. In the previous work with Mada and Tokihiro [8], the author 
presented a way to construct the conserved quantities of equation (12.41) which can also 
be applied to the negative solitons, by using certain gauge transformation to the BBSC. 
We described the structure of the conserved quantities for the BBSC, in terms of arclines 
connecting balls and vacant boxes. 



3.2. Coupled form of BBS(C). Hereafter we consider the BBSC with box capacity L 
and carrier capacity /. 



Lemma 3.1. Equation (12. 6 p is equivalent to the following coupled equations 

'Ul+^ = min(c*, L-Ul) + max(0, Ul + c\- I), 

with the boundary condition 



(3.1) 



lim c- = lim U- = 0. 

i—^—oo i— >— oo 

Proof 

If we determine c* by c* = EI=-oo(f^fc - ^k^^) in UM, we have □ 

On the contrary, we have c- = J2^k^-ooi^k ~ ^k^^) from the last equation in (13. ip with 
the boundary condition. Then the first equation gives (12. 6p . 

Remark 3.1. The equation (13. ip is also called an "ultradiscrete Yang-Baxter map" . 
Corollary 3.1. The coupled equations of the BBS are expressed as 



(3.2) 
Proof 




- min(c*, L - Uj) 



u + d-u 



We have only to set / = +oo in (13.11) . □ 
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n+1 



Cn+2 



t+1 



TTt + 1 



,t+l 



-n+1 



,t+l 
'n+2 



Figure 1. The time evolution of the coupled BBSC equation 

3.3. Periodic BBS(C). We consider the time evolution of the BBS(C) with periodic 
boundary condition. For a detailed discussion of the periodic BBS with Ui, Ci being non- 
negative, see |H] and [10]. We extend the results in P HO] to the periodic BBSC that 
allows negative solitons. 

Definition 3.1. We say that the set of initial values U\, f/g? ' ' '^n evolves as a periodic 
BBSC of size N, if and only if there exists at least one c\ such that c^+i = c\ and such 
that Ul'^^, 1/2'^^ , ■ ■ ■ U^j^^ are defined uniquely, independent of the choice of possible c\ 's. 
□ 

Remark 3.2. For some initial values {U^}fLi, there may exist more than one c\'s such 
that c^+i = c\ and such that {U^'^^} depend on c\ (Example \3.3]) . The time evolution of 
the BBSC is not unique in these cases, and we exclude such cases for simplicity in this 
paper. 

3.3.1. Periodic BBS. First, we deal with the periodic BBS f l3.2p . 
Lemma 3.2. 

c*;v+i = max[c\,vl] + N 

Here, 

v\ = max [iL - 2{Ul + ■ ■ ■ + Ul^) - U',], 
N = 2{U{ + U'2 + ---U%)-NL. 

( See figure [H ) 
Proof 

First note that N does not change under the time evolution of the periodic BBSC. We 
omit the superscript (■)* of cf s and f* for convenience and consider ctv+i =: /(ci) as a 
function of Ci. 

Fix one Ci such that Ci > Vc- 

Since ci> Vc> L — Ui, 

Ul+^ = min(ci, L-Ui) = L-Ui, 
C2 = ci + 2Ui - L. 
The inequality ci > Vc > 2L — 2Ui — U2 gives 
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f(Ci) 




Cl + N 



V, + N 



Figure 2. 



The graph of cn+i = /(ci) 



C3 = C2 + 2(t/i + f/2) - 2L. 
Using Vc> iL — 2{Ui + ■ ■ ■ + t/j-i) — Ui, repeated calculations lead to 



holds. Hence {Ui,Ci) does not depend on ci if i > J + 1. In fact {Ui,Ci) will be a 
constant that only depends on the initial values (f/i, ■ ■ • , Un)- Therefore, cn+i = /(ci) 
is also a constant if ci < Vc- It is easily seen from (13.21) that /(ci + 1) — /(ci) = or 1, 
and we already know that /(ci + 1) — /(ci) = 1 is equivalent to Ci > Vc- Thus we obtain 
/(ci) =Vc + Nifci< Vc. □ 

Theorem 3.1. The BBS (13.21) evolves as a periodic system if and only if 



holds for the initial values {U^y^i. 

Proof 

From the lemma 13.21 

• If iV > 0, no c° G Z satisfies /(c°) = c?. 

• If iV = 0, we have /(c?) = c? for all c? > Vc- 

• If < 0, the only c? that satisfies /(c?) = c? is c? = + iV. 



Q+i = Ci + 2(f/i + f/2 + ■ ■ ■ + f/i) - 



for i = 1, 2, ■ ■ ■ , A^. (This is easily seen by induction.) 

Thus we obtain c^+i = /(ci) = Ci + iV if ci > Vc- 

If Cl < Vc, on the other hand, there exists 1 < J < N such that 




N <0 
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In the second case, Ul is independent of the choice of c° because we know from the 
previous lemma that f// = 1 — (1 < < A^) for any > Vc- The evolutions are the 
same in t = 2, 3, ■ ■ ■ . □ 

Remark 3.3. Although normally we treat the case where Uf, c\ G this proposition is 
also valid for arbitrary real values ofU^, c\. □ 

3.3.2. Periodic BBSC. Next we deal with the periodic BBSC f l3.ip . As we are investigat- 
ing whether the system evolves from the time t to t + 1, we can omit the superscript (■)* 
of the variables below. 



Theorem 3.2. Let, 
M 
cl 



N. 



odd 



2(f/i + ■ ■ ■ + Un) - NL, 

max [L - 2{Ui + ■ ■ ■ + Ui) - a,+i + iL], 

i€{0,l,-N~l} 

min \l - 2(Ui H h f/i) - t/j+i + iL], 

ie{o,i,---N~i} 

{i\l<i< N, i: odd}, 
{i\ 1 < i < N, i : even}, 



d, 



df 



max 



mm 



max (/ — Ui), max {Ui — L + I + 1) 

even 



min {L-Ui- I] 



min {Ui 

even 



+ 1. 



When we put cat+i = /(ci), the following cases occur. 

• If oo > I > L: 

— If M ^0, there exists a unique Ci G Z such that /(ci) = Ci holds. 

— If M = and cl < cr, we have /(ci) = ci for all ci with < ci < cr. 

— If M = and cl > cr, there exists a unique c\ such that /(ci) = ci holds. 

• Ifl< L: 

— If N is an odd number: 

There exists a unique Ci G |Z such that /(ci) = Ci holds. 
(For Cl G |Z \ Z, we do not have a normal box ball interpretation. However, 
Vi Uf E Z ^ \Ji Ul^^ G Z holds, which means the number of balls in the box 
is an integer all the same.) ■ ■ ■ (**) 

— If N is an even number, 

* For di > dji, there is a unique Ci G Z that satisfy /(ci) = Ci. 

* For di < dji, we have /(ci) = Ci for all Ci such that d^ < Ci < dn, but 
{C/*"*"^} depends on ci. 

(This is the only case where we cannot determine the time evolution 
uniquely.) ■■■(***) 

• Ifl = L, there is a unique Ci that satisfies /(ci) = Ci. 

□ 
Proof 
If / > L : 

We regard Cj+i as a function of q to find 



+ 1) ~ Cj+l(Cj) 



1 {L-Ui<Ci<l-Ui] 
(otherwise), 
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which leads to 

'l (V«, L-Ui<c^<l-Ui), 
(otherwise) 

li L — Ui < Ci < I — Ui for all z's then, by induction we have 



/(ci + 1) - /(cO 



Ci+i = ci + 2(f/i + ■ ■ ■ t/i) - (i = 1, ■ ■ ■ , AT). 
Thus the condition (Vi, L — Ui < q < I — Ui) is equivalent to the following: 

(3.3) (Vi) L - 2{Ui + --- + Ui)-U+i+iL<ci<l- 2{Ui + ■ ■ ■ + U) - U+i + iL. 
The inequality (13.31) can be expressed as 

cl < ci and ci < cr, 

and we have 

[ /(cl) (ci < cl), 
cat+i = /(ci) = < Cl + M (cl < Cl and ci < cr), 

(ficR) {CR<C,). 

(Note that cl niay be larger than cr.) 

By examining the intersection of y = /(ci) and ?/ = Ci , we obtain the desired result. 
(When M = and cl < cr, possible ci's are not unique, but we have f/*"^^ = L — Uf 
regardless of the choice of Ci which makes the time evolution of f/* unique.) 

We conclude from the above argument that the system uniquely evolves in time ill > L. 

Iil<L : 

We have 

il-Ui<Ci< L-Ui), 
(otherwise) . 



Cj+l(Cj + 1) — Cj_(.i(Cj 

Hence, 



1 (Vz, I -Ui < d < L -Ui) and (A^ : even), 
/(ci + 1) - /(ci) = { -1 (V?, l-Ui<Ci<L-Ui) and (iV : odd), 

(otherwise). 
If I — Ui < Ci < L — Ui for all i, we have 

f/*+i = 2Q + f/*-/, 

for each i. Thus we have 

Cl {i G A^oddj 

/ - Cl (i e Nc, 

Therefore, the condition (Vi, I — U < Ci < L — Ui) is equivalent to the following set of 
inequalities: 
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Cn+1 



'N+1 



dL dR 
N:odd 



-ci 















dL dR 
N:even 



Figure 3. Graph of cn+i = /(ci) of BBSC under the condition / < L 



^ \Ui - L + I < Ci < Ui {te iVeven). 

Thus if N is an odd number we obtain 

f/K) (Cl<rfL), 
/(ci) = < / - Ci (di < ci and Ci < (Ir), 
(fidR) {dR<ci). 

(See the left graph of figure [31) 

Therefore, there is a unique ci such that /(ci) = ci (ci G |Z). Even though ci may 
not be an integer, as long as {f/*} are all integers, {U^^^} are closed in Z. In particular 
if di > dji, /(ci) is constant in Ci. 

If is an even number, we obtain: 

{/(di) (rfi < dL), 
ci {di < ci and ci < d^j), 

/(rf/j) (rf/j < Ci). 

(See the right graph of figure [31) 

Therefore if di > dji, /(ci) is constant in ci, which renders the time evolution unique. 
If di < dji, by fixing ci such that d^ < ci < d^ we obtain 

t/*+i = f/* + (-l)^(/-2ci). 

This indicates that the time evolution is dependent on ci. We do not have a unique 
periodic BBSC in this case. □ 

Corollary 3.2. The time evolution of the BBSC (13. ip with the periodic boundary condi- 
tion is not well-defined if and only if the three conditions 

I < L, 

the system size N is an even number, 
dh < dji, 
hold simultaneously. 

Example 3.2. We show an example of the periodic BBSC. We express the time evolution 

V = min(c, L - U) + max(0, U + c-l), 
d = U + c-V, 
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as follows: 
U 

c-l-d. 
V 



Let N 




8 


L 




1 


and I ^ 1 
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Example 3.3. We present two irregular cases in defining periodic BBSC. Here is an 
example of the case (**) in proposition \3.S\ Let N = 7,L = 5,l = l and b = j. 

14 14 14 1 
b- I -b- I -b- I -b- I -b- I -b- I -b- I -b 
14 14 14 1 

Next we show an example when the periodic BBSC is not well-defined. Let N = 4, L = 5 
and 1 = 2. 

2323 2323 
1-1 -1-1 -1-1 -1-1-1 , 2- I -0-1 -2- I -0-1 -2 
2323 4141 

We have more than one types of time evolutions depending on the choice of Ci . 

Remark 3.4. The BBSC is an invertible system, which is true even when some of the 
variables {f//, c*} take negative values. 

3.4. Periodic BBS with K kinds of balls. We consider the extended BBS where we 
have K kinds of balls distinguished by integers 1 < k < K. We call this system BBSi^. 
Here we consider the case where the capacity of the box is 1 in each box. Let U^^ be the 
number of balls k in the i-th box and c* the number of balls k in the carrier when the 
carrier is located between the [i — l)-th box and the i-th one. The time evolution rule of 
the BBSft: is 

(fc-l K \ 

j=i j=k J 

where k = 1,2, ■ ■ ■ ,K (See figure H]). 

Definition 3.2. The BBSk can be defined as a periodic system if there exists at least one 
set of 
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(U*i.Uu.-.U'k) 



('^f+l.l' '-f+1,2' '" ■ '^f+l.K) 



Figure 4. Time evolution of BBSi^ 

such that c^+ij = c* ^ hold for all j G {1, ■ ■ ■ , K} and the Ul\^ 's are determined uniquely 
independent of the choice of possible c\ 's. 

Theorem 3.3. The periodic BBSk is well-defined if and only if 

N / N K \ 

\ / 

holds for all / G {1, 2, ■ ■ ■ , A^}. 
Proof 

First note that we only have to consider the case where t = because Uj i is 

independent of t. We know from the time evolution rule that c^+i ^ depends only on 

1) 2) ■ ■ ■ ! c'j' ^ and does not depend on c\ ^^j^^, c\ ^,^35 ' ' ' ^^1 k- We omit the superscript 
(■)° for convenience hereafter. 

For / = 1, 2, • • • , iV let 

N K N 

j=i k=i j=i 

The condition in the proposition is rewritten as A*"; < for all /. 

We first determine the cases that allow for a ci^i such that satisfies c^+i^i = Ci^i. We 
calculate in the same way as in lemma 13.21 to obtain 

CAr+1,1 = max[ci,i,7i] + A^i, 

i-l i K 

j=i j=i k=i 

Thus there exists Ci^i such that C7v+i,i = Ci^i if and only if Ni < 0. 
Next we determine from ci_i and ci^2- We obtain 

Cn+1,2 = max[ci,2,72] + ^^2) 

N / K \ N 

j=l \ k=2 J j=l 

and 72 is a constant determined by the Ui^kS only. Therefore we have ci^2 such that 
CAr+1,2 = ci,2 if and only if iVg < 0. We know from ^1,1=1 = J2i=i that N2 = N2. 
By repeating the same calculations for CN+i,k {k > 3) we obtain the proposition. If Nj = 
for some j there exist more than one cij such that cn+ij = Cij. The time evolution of 
Ui^k is still unique in this case. □ 



7i = max 

ie{i,-,N} 
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3.5. Some elaborations on conserved quantities. 



Proposition 3.1. The invertihle transformation 




{(Jl^^ = min(5*, (L + 2m) 



Ul) + max(0, Ul 



+ c\-{l + 2m)) 



Here m is an arbitrary real number. 



Proof 



We immediately obtain fl3.4|) by direct calculation [TT]. □ 

When m > 0, in particular, we call this transformation an "upward-shift translation" . 
If some of the variables take negative values in the initial conditions, by putting m > 
large enough, all ?7*'s and c-'s are transformed to be positive. Note that the capacity 
of the boxes and that of the carrier increase by 2m, and that the boundary conditions 
will change: Uf is not at i — i- ±oo anymore, but it becomes some nonzero constant at 
i — )■ ±oo. Also note that the value [2(?7° + [/^ + ' " " U%) — NL] does not change through 
this shift. Hence, treating the negative solitons is equivalent to treating the non-negative 
BBS's with the boundary conditions that the solutions [/^ have some positive constant 
value at 72 — )■ ±oo. 

Remark 3.5. In this transformation, the values Ul at i ^ +oo and those at i ^ — oo 
have to be the same. In section 5, we present a way to deal with the boundary condition 
such that limi_>_oo f/* 7^ limj^+00 . □ 

The construction of conserved quantities found in [8] can also be performed for a peri- 
odic BBSC with general L and /. We depict by arclines the movement of the balls from 
the box to the other box according to the time evolution of BBSC. The operation of 
drawing arclines when the capacity of the carrier is / is denoted by (OP)/. 



• We take out balls from the boxes from the left to the right according to the 
evolution of BBSC. We distinguish each ball and when two or more balls are 
taken out from the same box at each step, the one at a lower position is taken out 
first. The carrier can hold at most I balls. 

• When two or more balls are passed from the carrier to a box at the same step, the 
one taken by the carrier last will be the first to go back to a box. 

• We depict by arclines the movement of the balls from the box to the carrier to the 
other box. 

We can see some structures in the set of arclines connecting the balls and the empty 



Theorem 3.4 (Kanki-Mada-Tokihiro [8j). On the periodic BBSC upward-shifted by m, 



time evolution of the system. We also denote by ((7°, , ■ ■ ■ ) the set of these constants 



(OP)/ 



boxes. 
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Figure 5. Drawing arclines to the periodic system of size = 12. 

for the "vacuum state" where every box has exactly m balls. The difference of these two 
{Ci — C°, C2 — C2 ) ■ ■ ■ ) ^■^ ^^■^^ '5^^ conserved quantities of the BBSC. We rewrite it 
(Ci, C2, ■ ■ ■ ) and call it "the set of conserved quantities of the BBSC" . □ 

We extend this theorem to the general periodic BBSC. 

Proposition 3.2. For the periodic BBSC with parameters L and I and initial conditions 
{[/°} that satisfy the following 

For both the initial conditions {U^} and the vacuum solution {m} 

we can define the periodic system for I G [L, +00], 

(3.5) and also for all I, there exists Ci G Z such that /(ci) = Ci, 

we can apply theorem \3.4\ to construct the conserved quantities. □ 

In other words, we consider the cases other than (**) in proposition 13.21 Note that for 
I with I < L we allow the case (* * *) in proposition 13.21 

Example 3.4. See figure O for an example. We consider the periodic BBSC with a 
box capacity 1 and system size N = 12 and set the initial value to be 001101101000. 
By an upward-shift with m = 1, initial values are transformed into 112012212111. We 
have iC^,C2,C,,C,) = (6,5,3,1), (C?, C^, C3O, QO) = (6,6,0,0) and (^1,^2,^3,^4) = 
(0,-1,3,1). 

3.6. Defining Background Solutions. We propose the following way to distinguish the 
background solutions from the positive solitons. Note that the balls connected in (OP); 
are included in those connected in (OP);/ where /' > I. 
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Figure 6. (Example 13 .Sp : Distinguishing between -11-1 and -12-1. 

Definition 3.3. For the BBS upward- shifted by m, the balls connected in the process 
{0P)2m constitute negative soliton solutions. Remaining balls constitute positive solitons. 

Definition 3.4. The solution {f/*} is a "background solution" if Ck = for^k > 2m + 1. 

Example 3.5. For the BBS with the box capacity one, ■ ■ ■ 01110 • • • and ■ ■ ■ 01210 ■ ■ ■ are 

both stationary solitary waves moving at speed one. We can distinguish these two waves 
by constructing conserved quantities. The former is a background solution without positive 
solitons, and the latter can be interpreted as a background solution with a soliton of length 
one superimposed. We shift both systems upward by m = 1, and find that the former has 
the conserved quantity {Ci,C2) = (0, —1), Ck = 0{k > 3). On the other hand the latter 
has {Ci, C2, C3) = (0, -1, 1). (See figure\B\) 



4. Relation to the discrete system 
We consider the following transformation 

(4.1) cr* = r(*-")V2V* 

for the solution o"^ of the bilinear discrete KdV equation. 

Proposition 4.1. The transformed function satisfies the following gauge transformed 
discrete KdV equation 

(4.2) (1 + SKWal,-' = 5^'^''^^~al_,\~<^' + ~a'^aU,. □ 

This statement is proved by a direct calculation. We consider the ultradiscrete limit 
of the equation (14. 2p . We let 5 = e~^/^, a^^ = e"^"/^ and = e'^"/^, and the following 
proposition holds. 

Proposition 4.2. For f/* and f/* defined as U^^ = + t^Ii — T^^^ — t^_i and f/* = 
+ ^n-i ~ ^n^^ ~ ^n-i? relation (14.11) is transformed into 

=U^ + 1. 

n n ' 

in the limit e — )■ +0. □ 
Proof 

We take the limit e — -|-0 to obtain f* = r* -|- |(t — n)^. Therefore we have f/* = f/* + 1. 

□ 
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Remark 4.1. The transformation (14. ip of the discrete KdV equation corresponds to the 
'1 ' upward- shift to the BBSC in the ultradiscrete limit. Thus the solution of the equation 
(14. 2 p goes to a '1' upward- shifted solution of the BBSC in the ultradiscrete limit. 



5. Background solutions from Jacob: theta functions 
We define Jacobi theta function by 

oo 

n=— oo 

wliere q = e™'^, z = e^'^'" and i = a/— 1. We suppose lm{ri) > 0. Tlie Jacobi tlieta 
functions are quasi-doubly periodic functions, that is, we have 

Proposition 5.1. Let v = ri{ni — n2 — ^3) + r/o- Then the function T{ni,n2, n^) '■= '&3{v) 
is a solution of the following equation, which is a gauge-transformed discrete KP equation 

(5.1) (sa)rir23 - r2ri3 + (1 - s)nTi2 = 0. 

Here, a = e^^'"^{= g^) and s is an arbitrary complex number. The lower script 'i' 
of the function r denotes a '+1' shift in the variable ni, i.e., ti = r(ni + I,n2,n3), 
T23 = T{ni,n2 + 1, n3 + 1), etc ■ ■ ■ . 

Proof 

From the quasi-doubly periodicity we have 

n = ^3(t; + r7) = e-'"(2.+^)7?3(^), 

r23 = Mv-'^v) = e'''^'''"-'"^Mv), 

which lead to 

(5.2) r^T2s = e'''^'^-'^\Mv)Y. 
We also obtain 

(5.3) r2ri3 = e^'^(2^-'')(^93(t^))', 

(5.4) r3ri2 = e*"(2.-,)(^^(^))2_ 

Thus we have the result. □ 

Remark 5.1. // a function a{ni,n2,ns) satisfies the relation aTiT23 = T2T13 = T3T12, 
then r(ni,n2,n3) satisfies the same equation (15.11) . The discussions in this section are 
therefore equally valid for such functions a. 
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5.1. Reduction. Since the parameter s is arbitrary, we can put s to be s = e*'^^-^ > 0). 
The parameter C, is related to the box capacity of the BBS. From the way v is chosen, 
we have r = ri2, which is the reduction condition from the KP equation to the KdV 
equation. If we rewrite := r(t, 0, n) the following lemma holds. 

Lemma 5.1. The gauge-transformed discrete KP equation f lS.ip is reduced to the following 
gauge-transformed discrete KdV equation 



(5.5) SarJ' '^n+l + (1 ^)'^n+l'^n — '^n ''"n+1- 

Remark 5.2. We have the following relation with the upward shifted discrete KdV equa- 
tion 



1 + 5' 
a = □ 

Although the coefficients in equation (15. 5p are different from the bilinear form of the 
normal discrete KdV equation, we find from the following proposition that this does not 
lead to the loss of generality. 

Proposition 5.2 (Tsujimoto-Hirota ^). By change of variable 

t-i 

^ 'n'n+l 



'n+1 'n 



the equation (15. 5 P is transformed into the discrete KdV equation 

"^n+l 

where 6' = sa. 
Proof 

From the following identity 

(LHS of dEDn^n+i XT* - (LHS of dESD) X r*+2 
= (RHS of m))n^n+i XT* - (RHS of (153])) X r*+2, 



we have 
Multiplying 



^t+l t t-l _ t+l t t-l _ r/ t t+l i-l _ r/ t t-1 t+1 
^n+2'n'n+l 'n+1 'n+2'n ~ " '^«'n+l'n+2 " 'n+2'n+l'n ■ 



'n+1 



on both sides, we have the result. □ 



'^n+l''^n+2''n''n+l 



5.2. Coexistence of Solitons and Jacobi theta backgrounds. If we ultradiscretize 
the solutions of the equation (15. 5p . we will obtain the solutions of the BBS. We will 
ultradiscretize the equation and the solutions. To do this we first assume that r] and rjo 
are both purely imaginary and then replace inr] with —L/e. (This choice of e is empirical.) 
Then we take logarithms on both sides and take the limit e — )■ +0. 

Note that our paper is not the first to ultradiscretize the theta functions. For example 
in [12] , Iwao and Tokihiro ultradiscretized the theta function solution of the periodic Toda 
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equation. Our method is unique in that we connected the ultradiscrete theta function to 
the background solution of the BBS. 

Lemma 5.2. The ultradiscretization of the solution = i^slv) is the tau function of 
the constant solution which takes = 2L for n e {—oo, oo). Here we have the following 
relation between Ul and 91 . 



Proof 



exp [—iTi7]k{2{t — n) — k + 2rio/ri}] . 



From the first hne to the second hne we replaced k with —k. We then transform iirr] into 
— J, and affect on both sides lime_j,+o e log(-). If we write the left hand side 9^ after the 
ultradiscretization, we have 

ei = m^x[Lk{2{t-n)-k + r]')]. 

Here we put rj' = 2?7o/'7- The solution constructed from 0^ takes constant value 2L 
on n e {—oo, oo). □ 

We now have the background state '&3{v). We can add A'^-soliton solutions onto the 
background state ^3{v). If we suppose that r* := -(93(^)0^ is also a solution of (15. 5p . then 
0^ satisfies the ordinary discrete KdV equation 

s0^+V*:;\ + (1 - = 0n'V^i\- 

We can take 0^ to be the A''-soliton solution of the discrete KdV equation. We denote 
the ultradiscrete limit of 0*„ by $^ and define Vj; by V„* = ^li + ^t^-i - ^rt^ - ^n-i- Since 
s = e^'^^^ = e~^^/^, the solution Vj^ corresponds to the solution of BBS with box capacity 

Proposition 5.3. We have the BBS with box capacity (4 + C)L from the equation (15.51) 
through ultradiscretization. 

Proof 

The proof is based on Takahashi and Hirota [TS]. If we put 6' = se^^'"^, then (15.51) is 
transformed into 



from the proposition 15.21 Thus we have 



From (15.61) we obtain for any M < 
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Here we have 

t _ 'n'n+l _ ^n^n+\ VnVn+1 

'^n — t t-1 ~ M Qt-1 Mt At-1 ' 
'n+l'n '^n+l'^n S^n+lVn 

where 

and the positive sohtons satisfy hm„^._oo0^ = 1- 

Thus we obtain hm„_!._ooW^ = e"^*'^''. Therefore in the hmit M — +oo we have 

" w'+\ 1 



k=—oo 



which is equivalent to 



w. 



k=—oo ^ 



After replacing mrj with — j, and with e^"/^, we take the limit £ — t- to obtain 

n 

UiX\ = - max[[/*^i - (4 + 0^, E (f/f' - Ul)]. 

From — max(a, h) = min(— a, —b), the equation is the BBS with the box capacity (4 + 
OL. □ 

These results can be summed up to the following theorem. 

Theorem 5.1. As solutions of the gauge-transformed discrete KdV equation (15.51) . we 
have the following type of solutions constructed from the Jacobi theta functions. 

{Background solution) {N-soliton solution) 

The ultradiscretization of this solution is a solution of the BBS with box capacity (4+^)L 
and IS 2L + 1/* . □ 

Example 5.1. For example when L = 1/2,^ = 2, we have a solution of the BBS with 
box capacity 3 as shown in figure \2\ 

We successfully obtained the positive solitons superimposed over the background 
solutions U^. However, to obtain negative solitons we have to consider solutions with 
holes in the background solutions, which needs a more detailed discussion [S]. 

6. BBS WITH IRREGULAR BOUNDARY CONDITIONS 

Finally we present the BBS with a boundary condition where 

lim f/*7^ lim f/*. 

n— >+oo n— >— oo 

Lemma 6.1. The t- function of the N-soliton solution of the BBS is as follows 

(6.1) T*=max [j;P,jt-(^min(L,P,)]n+[^0i]- J] 2mm(P„P,0 
Vie,/ / VjgJ / \ieJ / i>j ijeJ 

where Vi Pj > and Li G M. 
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Figure 7. Example of background and i-soliton solution constructed from 
Jacobi tlieta function 

We take L = 3 and consider the soliton solution made up of solitons {Pi,Oi) = (2,1) (z = 
1, 2, ■ ■ ■ , A^o)- It has the expression 

t(x)= max [2|J|a; + IJI -2|J|(|J| - 1)1. (x = t-n) 

0<|J|<Aro 

In the limit Nq oo, the A^^o-soliton solution converges to the following background 
solution: 



(6.2) 



Ugix) 



max 

fc>0 



2k [x + -] ~2k{k- 1) 



fo 



Values of f/* = Ug^x) + Ug^x + 2) — 2Tbg{x + 1) at integer points at time t = are as 
lows. 



n 


■■■0 12 


3 4 5 ■■■ 


^ n 


1111 






The 



Dackground so^ 



ution (16. 2p has the new boundary condition where the solution tends 



to 1 in n —7- — oo and to in n — )■ oo. We can add positive solitons to (16. 2p . 

Example 6.1. We construct the following 2-soliton solution travelling in the background 
state (16. 2p . Let the width and the phases of the two solitons to he added he 



Pi = 4, ^1 = -8; P2 = 5, 02 



-12, 



and we denote these solitons hy soliton A and B respectively. The general form of the 
solution t{x) is given as 

r(x) = ma.x[Tbg{x) , 4t — 3n — 8 + Tbg{x — 2), 

5t - 3n - 12 + Tbg{x - 2), 9t-6n-28 + ng{x - 4)]. 

Construction of the coexisting state of solitons and negative solitons is hased on [8J . What 
happens during the time evolution of this system is as follows. (See also the time evolution 
pattern of the system at the end of this section.) 

• < t < 5 
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We observe that the soliton A travels at speed 2 and the soliton B travels at 
speed 3 on the background 1. The soliton B takes over the soliton A in the same 
way as in BBS with box capacity 1. 

• i ~ 6 

The soliton B climbs down the discontinuity of the background state from 1 to 
0. The phase of the background is shifted 2 to the left. The soliton B travels at 
speed 5/3 from now on. 

• i ~ 12 

The same incident happens to the soliton A and it travels at speed 4/3 hereafter. 
The phase shift of the background is also —2. 

Remark 6.1. The BBS may take different boundary values at n ^ oo and n — )■ — oo 
like in this example, which is not treated as a periodic BBS. Other irregular boundary 
conditions can also be considered and general solutions are calculated in the same way. 
In this section we have only dealt with the case where L — 3. We can consider general L 
to obtain the boundary condition in which the left side and the right side take arbitrary 
integer values. 



t= 
t= 
t= 
t= 
t= 
t= 
t= 
t= 

t= 8 

t= 9 




1 
2 
3 
4 
5 
6 
7 



:222: : : : 
: : :222: : 
: : : : :222 



22 : : : : : : 
:22: : : : : 
: : 22 : : : : 
222: :22: : : 
: :22: :222: 



:::::::: :22: : :222: :. . 
:::::::::: 22 :::: 222 . . 
:::::::::: :22: : : :131. 
:::::::::::: 22 : 
:::::::::::: :22 
t=10 ::::::::::::::: :22 
t=ll :::::::::::::::: :22 
t=12 ::::::::::::::::: :3 



.23 

. .32. . . 
. .131. . 
. . .23. . 
32. 



t=13 ::::::::::::::::: .31. . . 

t=14 ::::::::::::::::: .22. . . 

t=15 ::::::::::::::::: .13. . . 

t=16 :::::::::::::::::. .31. . 

t=17 :::::::::::::::::. .22. . 

t=18 :::::::::::::::::. .13. . 

t=19 :::::::::::::::::.. .31. 

t=20 :::::::::::::::::.. .22. 
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.23 

. .32 

. .131. . . 
. . . 23 . . . 

32 . . 

131. 

23. 



The time evolution of 2-solitons on an irregular background is shown above. Here the 
frame itself moves along with the background state so that the discontinuity in the back- 
ground solution seems fixed. Note that the background state itself is moving to the right 
at speed 1. The symbol ':' indicates a background 'i' and '.' a zero respectively. We see 
that the phase shift of the background after colliding each soliton is —2. 



7. Concluding Remarks 



We first discussed the conditions under which the periodic BBSC with general box and 
carrier capacities is well-defined. We extended the construction of the conserved quantities 
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to general BBSC. Conserved quantities are useful in distinguishing background solutions 
from positive soliton solutions. We then showed that the Jacobi theta function is a solution 
of the gauge transformed discrete KdV equation, and that the ultradiscretization of this 
solution corresponds to the background solution of the BBS. The author wishes to extend 
this method to more general functions in order to deal with negative solitons. Finally, the 
BBS with irregular boundary conditions has been constructed. To obtain the solution of 
the system with general L is a future problem. 
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